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HOMMCLATIIRE 

A1 * Top diameter of the tapered rod, inches 
Cl Bottom diameter of the tapered rod, inches 
ID Internal diameter of the glass tube 

D1 Internal diameter of the float, inches 

D2 Orifice diameter of the float, inches 

D3 Diameter of the tapered rod at a particular height of float, 
inches 

HI Height of the float, inches 

H2 Length of the tapered rod, cms. 

H3 Orifice he.ight of the float, inches 
H4 Length of the pipe support, cms. 

•HT 1 Bloat reading at a particular height, cms, 

Q Volumetric flow rate in cc/sec 

R Drag force 

Density of the float, gm/cc 
f Density of fluid, gm/cc 

Density of water, gm/cc 

P.j -Pressure at float-top 

P 2 Pressure at float Bottom 

Ap Pressure drop due to float 

&P* Pressure drop in the flowing liquid 

U Viscosity of water, gm/cc-sec. 

VS Volume of float, r% 

WS Mass of float, gms 
0 Taper angle of ro^ ^ egress 


* *. « 
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ABSTRACT 

An inexpensive and s - * reply constructed visual flowmeter 
has been designed. It consists of a. constant cross-section; 
glass tube in conjunction with a tapered rod and a float. 

The variable area is provided by the annular region between 
the tapered rod and the glass tube. The liquid flows through 
the float. Thus, this flowmeter is also capable of measuring 
flow rates of opaaue and translucent liquids. It has a 
linear calibration curve. 

Following experiments with water formulae relating 
the variables of the meter in the form of dimensionless groups 
have been developed. These can be used to specify the dimensions 
of the motor for a given duty. 



CHAPTER I 


INTRODUCTION 

Tlie measurement of incompressible fluid-flow is often 
made with the help of variable-area flowmeters. Of these, - 
the rotameter is the most commonly used in industry. The 
precision needed in industry is quite adequately covered by 
a rotameter? besides, its ether advantages are ease of 
handling, direct readability and an easy access for cleaning 
purposes. 

Presently in India, rotameters have a high initial cost 
and even subsequent replacement of the glass tapered-tube is 
not cheap. This together with the difficulties encountered 
in reproducing the taper on the glass tubes, and hence dupli- 
cation of calibration , limits the desirability of this type 
of flowmeter. 

The aim of this study is to explore the possibility of 
replacing the tapered tube with one having a constant cross- 
section. A concentric tapered rod fixed at the centre of this 
tube would provide the variable area for the flow measurement. 

A float; concentric with this central tapered -rod would slide 
along the inside of the glass tube. It is estimated that such 
a flowmeter can be manufactured with ease and at low cost 
because of the fact that the highly sophisticated and expensive 
process of manufacturing tapered glass tubes has been eliminated. 
Moreover , re placement costs of constant cross-section tubes for 
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this flowmeter will be nominal, and duplication of calibration 
also simple. 

Experimental observations were made in order to establish 
expressions relating the parameters that affect the performance 
of the flowmeter, and to test the feasibility of operation.- The 
method of dimensional analysis was chosen, from among the 
various methods, because it helps in producing dimensionless 
groups, which are in themselves lesser in number than the 
existing variables. Furthermore, it also enables to obtain 
dynamically similar systems, which are needed for scale-up 
processes. This method does not in itself provide a complete 
solution, rather it yields functional relations which become 
explicit after combining with appropriate experimental data. 

These equations, can then explain the effect of various quantities 
upon the performance of the meter. 





CHAPTER II 


LITERATURE REVIEW 

The variable area flowmeter was first invented by 
Edmund Augustine Chemeroy of Paris in 1868. He used a tapered 
pipe and a float having mechanical transmission to read the 
flow-rate. Later on, it was repatented several times in 
different countries between the years 1880 and 1905. 

In England, Beacon got a patent in 1875 for a cone-and- 
disk type flowmeter. Then Alfred Ewing in 1879 for the first 
time, used a tapered glass tube for measuring liquid flow-rates. 
George Jos tin, obtained in the same year a patent for measuring 
gas flow rates. 

Clausolles, of Prance, in 1903, found a velocity meter 
designed to measure velocity of a liquid flowing in a pipe, 
using a tapered tube and a float. 

Karl Kupper was the first to use the word rotameter for 
variable area flowmeter in 1908 because of the rotary motion 
of the float resulting from the inclined slots made on its body. 
Since then this name has been retained for this instrument 
inspite of the various changes it has undergone till now. 

In 1934, Bentzel used a Darcy type tube, one part of 
which was tapered, and the float indicated the water velocity 
in the stream according to the rating of the instrument. An 
excellent review of the historical development of variable 




5 


area flowmeters during 1868-1934 is given by Kolupailal 1 ^ 

( 2 ) 

Stout and Rowe constructed a flowmeter in 1938, 

consisting of an ordinary glass tube, a metal cone and an 

annular bob. They studied the effect of shape of the bobs 

shown in Figure 1, specific gravity of bob, taper of the cone 

and the changes in the annular area between the bob and the 

cone. Nearly in all the cases mathematical relationships 

were not obtained. 

(■*>) 

Fischer ' in 1940, designed a float which eliminated 
the effect of small changes in viscosity on the rotameter 
calibration. Hence he called it a stable-vis float. This 
f loai|consisted of a disc having a sharp-edge orifice, which 
served as a means of flow restriction. A separate body located 
in a pocket of stagnant liquid,’ attached to the disc through 
a tube of small diameter, helped in giving the necessary weight 
to it. This body , completely removed from the fluid flow? stream, 
did not provide extra drag. A central guide wire was used to 
guide this float. 

Kitner^ in 1942, developed a visual flowmeter having 
straight walled tubes. The glass tube had a concentric inner 
brass pipe which also constituted the inlet and outlet of the 
flowmeter. Two parallel-sided slots were cut along the 180° 
axis of the brass pipe. Above but connected with the slots, 
a number of holes were, cut through the pipe walls to allow 7 the 
fluid to return to the pipe after passing around the float. 
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The outer glass tube helped in confining the flow as well as 
enabled to read the float reading with reference to a scale 
inscribed on the slotted brt.s» pipe. Parabolic calibration 
profiles were obtained for different cases. 

/ r \ 

Coleman'''' showed that by making float density twice the 
mean fluid density, compensation for changes in the fluid 

density could be obtained. 

(6^ 

Danckwerts ' designed another constant cross-section 
variable area flowmeter in 1961. He used two concentric trans- 
parent tubes, the inner one being made of perspex and perforated 
by a regular array of holes. A rotameter-type bob was fitted 
closely inside the perforated tube, as a result of which most 
of the fluid passed out of the holes below it. The minimum 
flow which could be measured was that which leaked past the bob. 
The height of the bob indicated the flow-rate. The biggest 
disadvantage of this flowmeter was that the flow-rate did not 

have a linear relationship with float height. 

(7) 

Kehat ' extended the work of Danckwerts in 1964, to get 
a linear calibration curve. He replaced the perforated inner 
tube of Danckwerts with a brass tube having two parallel longi- 
tudinal slits, in which the float was allowed to move. The 
upper end of the inner tube was kept open as that in Danckwerts 
case. He was successful in getting a linear calibration curve 
for a small range only. 

(q) 

Kehat 's flowmeter was further modified by Kuloor v ; to 
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get a still better performance. The only modification he did, 
was to close the outer end the inner tube by a cap. Thus, 
the fluid was permitted to flow, only through the slit. The 
results obtained were encouraging. Linear relationships upto 
10 cm. float height for all combinations of liquids and float 
weights were obtained whereas Kehat’s flowmeter for the same 
dimensions would have given linear readings only upto 2 cm. 
of float height. 

In 1965, Tarish/''* produced a highly sensitive rotameter 
by using a convergent-divergent tube with two floats, one in 
each tube joined together through a thin rod. 





CHAPTER III 


EQUlJb?ISHT DESIGN AND EXPERIMENTAL PROCEDURE 

3.1. Equipment Design : 

The flowmeter consists of three main components, a glass 
tube, a concentric inner rod. and an annular float. A circular 
metallic plate was attached to both ends of the rod which 
served the purpose of keeping it. concentric with •- 
the glass tube. This in turn was made to rest on a. support 
made out of a thin pipe. Por the present work, rubber corks, 
having inlet and outlet pipes fitted in the centre were fixed 
at both ends of the glass tube. These also helped in keeping 
the assembly in position. Figure 2 shows the details of this 
flowmeter. 

Class Tube s A glass tube having a constant internal diameter 
of 0.985 inches was employed throughout the experiments. Its 
length had to be varied according to the total height of the 
internal assembly of the flowmeter. 

Tapered Rod ; Several tapered rods shown in Pig. 3a, made out 
of aluminium, were used in this work. The length (H2.), top and 
bottom diameters, (A1 and Cl) and the taper angle (9) were 
varied to study their effects on the performance ofjthe meter. 
Details of the dimensions of the rods used are given in 
Table I. 
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TABLE 

Is DETAILS 

OP TAPERED 

RODS 

Rod Ho. 

Approx. 

9. 

Cl 

A1 

Hj: 




■' '^ch. 

inch. 

Cel. 

1 

1 C 


0.456 

0.250 

15.0 

2 

0.75" 


0.402 

0.2 50 

15.0 

3 

0.50° 


0.349 

0.250 

15.0 

4 

1.00° 


0.590 

0.380 

15.05 

5 

0.75" 


0.510 

0.351 

15.0 

6 

1 .00° 


0.590 

0.250 

24.65 

7 

0,75° 


6.515 

0.230 

24.65 


Float s The float, though small in size, is an important compo- 
nent of this flowmeter. Properties like stability of the float, 
which have a direct impact on the accuracy of the flowmeter 
were considered for design purposes. Unaided floats try to 
rotate as well as tend, not to remain concentric. Hence, to 
provide mechanical aid for stability, the outer diameter of 
the float was made of a size just sufficient to slide inside 
the glass tube. To reduce the friction between the glass tube 
and the float, the area of contact between the two was restricted 
to both ends of the float as shown in figure 3b. Small vertical 
slots were made on this area to allow air, entrapped between the 
float: and glas§ tube, to escape. when the flowmeter is empty. 

In accordance to Heads work, a sharp-edged orifice 

was used to provide the necessary constriction to fluid flow in the 

( 2 ) 

glass- tube, following the conclusions of Stout and Rowe , this 
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orifict was aori: .sod at the bottom of the float in order to 
keep the centre of gravity as low as possible. 

The dimensions of various portions of the float were 
selected usm/: the following guide lines, based on the work 
of Hoad^^ for ro tame ter- type floats. 

3Z -= 0.831 
0.531 ^ HI 3i 
K3 = 0.01 31 + 105 

where DZ = orifice diameter 

3i = inner diameter of float 
El = height of float 
TI3 = orifice thickness 

The floats were made from aluminium, nylon, and teflon whose 
dimensions arc given in Appendix 1. 

Metallic Base Pla t o? Two circular metallic base plates made 
of brass, were screwed on t© the tapered rod, one at each end, 
to keep it concentric with the glass tube. The thickness of 
this plate was 0.25 inches. To minimise pressure losses of 
the flowing liouid, the base was made of the shape as shown 
in figure 3c. 

Pine Support ; To support the tapered rod .without choking the 
flow, a piece of thin metallic pipe having a diameter just 
sufficient to slide into the glass tube was employed at both 
the ends . 
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PIG-USE 3bs {EIIE PLOAi FIGURE 3c ; TEE METALLIC BASE 
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This pipe also helped in reducing the expansion effects, 
to he discussed later in detail. The different lengths of 
these pipe supports that v.ere used are 1.1, 3.6, and 4.7 cm. 
respectively. 

3.2 Experimental Procedure s 

Water was used as fluid to check the performance of the 
flowmeter. It was made to flow from an overhead tank to 
the flowmeter. The upstream of the flowmeter consisted of a 
hy-pass line. Two valves were employed to control and regulate 
the flow rate. Metal piping of l/4 !l O.D. was used throughout, 
for making the connections. A schematic flow-diagram of the 
above system is shown in Figure 4. 

To minimize fluctuations in flow, the use of reducers was 

fill 

avoided in the upstream section . ' A manometer connected 
across the flowmeter, and a pre calibrated standard rotameter 
in the downstream section served the purpose of indicating 
the total pressure drop and the flow-rate, respectively. The 
out-going liquid was recycled to the overhead tank. , With the 
available head a maximum flow rate of 3.5 l/min. was obtained. 
The range of flow covered was from 0.4 to 3.5 l/min. . The 
overhead tank was of a very large capacity which ensured a 
constant head over any single observation. Throughout the 
experiments the temperature of water was 26 + 2°C. 

About 40 different sets of observations were taken using 
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ricrns combinations of different float sizes, taper rods etc. 
each set, on the average, 13 readings were taken for 
Iferent flow rates. Repetition of flow measurements showed 
)d reproducibility. 



CHAPTER IV 


RESULTS 


The changes in flowmeter properties such as? sensitivity, 
stability, rotation etc., which were observed during experimentatio] 
have been discussed under qualitative results* Results relating 
the effect of variables quantitatively, through mathematical 
models, have been grouped in quantitative results. Dimensional 
analysis and multiple linear regression analysis which were 
employed are given in Appendix 4 and 5 respectively. 

5.1 Q ualitative Results ; 

a) Calibr at ion Curve s Linear relationships were obtained 
throughout the span between float height and liquid flow rate 
except, for a small initial portion ranging between 3 cm. to 
4 cm. of float height. Log-log behavior was evident in this 
range. As an illustration the calibration curve for Rod Uo.1 
and. Float Ho.1 is shown in Figure 5. 

b) Effect of Inlet and Outlet Pipe Diameters ; The total 
pressure drop across the flov/mcter was observed to increase 
for smaller inlet and outlet diameters. The calibration curve 
of the meter did not change in these cases. An illustration 
is shown in Figure 5. Pipes having 0.8 cm I.D. were used 
finally throughout the later experiments to match the other 
fittings of the set-up. The results are given in Table II. 
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1.2 1 .4 1 .6 1 .3 2.0 2.2 2.4 2.6 2.8 3.0 


Flow Rate in Ldters/minute 

FIG-URE 5 s 


. CFLIERAT10H FOE FLOAT 1 ROD 1 
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TABLE II: PRESSURE DROP FOR DlFE'RPF.Tvrn 1 TKrra'PS aatd 
OUTLETS' 

Plow- rate at winch the data. was taken — 3.04 litres/ minute 


Inlet end outlet 
pipe diameter in cm, 


Difference in Height, 
cm. mercury 


0.8 ■ 1.15 

c) Effect of Inner Assembly on the Pressure Drop Across the 
Flowmeter : Experiments showed that the major contribution 

towards the total pressure drop across the flowmeter was only 

due to the inlet and outlet pipes. The inner assembly had a 

nominal effect on pressure drop. 

The results are shown in Table III below: 

TABLE III: PEESSUBE PROP DUE TO DIFFERENT COMPONENTS 
OP THE METER ' ' “ ' 


Inlet and outlet pipe diameter 
Observation flow rate 
AH, for glass tube, only 


0.8 cm, 


3.04 litres/minute 
1.15 cm. mercury 


AH, for glass tube + metallic bases = 1.15 cm. mercury 


AH, for glass tube + metallic bases + 
tapered rod 


1.2 5 cm. mercury 


AH is difference in height of manometer reading. 

d) Effect of Height of Metallic Pipe Support: The height 


of the metallic pipe support had an effect on the linearization 
of the lower portion of the calibration curve. An optimum height 
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Initial reading of float = 0.0 era-. 



CALIBRATI01T FOR FLOAT 1 ROD 1 
FIGURE 6 s EFFECT OF PIPE IEIIGTH 
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of this support, with which an entirely linear calibration 
could be obtained, was found to exist. This height was 
different for different cases of the inner assembly and was 
found to depend on the taper angle (Q), the bottom diameter 
(Cl) of "the tapered rod and the slope of the upper portion of 
the calibration curve. A mathematical relationship of the above 
variables that was obtained is given in section 5.2 of this 
chapter. Figure 6 depicts the above effect as an illustration 
for the case of Rod No.1, Float No.1. The rest of the data 
are given in Appendix 2. 

e) Effect of the Bottom Diameter of the Tapered Rod ; For 
the same float and nearly the same taper angle of the rod, the 
slope of the calibration curve increased with an decrease in 
bottom diameter (Cl) of the rod. As already discussed in the 
previous section, it governed the determination of the optimum 
height of the pipe support. The flowmeter showed better 
performance with larger bdt'fccm diameters. 

f) Effect of Orifice Annular Opening s Observations showed 
that there was a particular limit to annular opening beyond which 
the flowmeter did not work satisfactorily. For (E2-D3) ^0.3, 
and Rod Nos. 2 and 3, light floats attained equilibrium after 

a very long time. A considerable amount of rotation was also 
observed in them. The above effect was seen to increase with 
a decrease in taper angle of the rod. Heavier floats did not 
however, show rotation and unstability to that extent. 
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g) Effect of Inner Diameter (Si) of the Floats floats of 
nearly the same weight, hut with smaller (El), showed rotation a 
little "before the limit of f,E2-D3 ) ^ 0.3 was attained, floats 
having El = 0.78 inches worked quite well. 

h) Effect of Taper Angle of Rod s; The lesser the taper 
angle of the roe the more sensitive was the flowmeter, figure 8 
shows the calibration curms of the cases using Rods Eos 1,2, 
and 3, with float i! o . 1 . Observations with other rods and floats 
arc given in the Appendix 3 , 

i) Top Taper on float ; There was no significant effect of 
the top taper of the float on the performance of the flowmeter. 
The taper was provided to avoid the sudden expansion of the fluid 
as it crossed the float. 

3) Pressure Erop Across floats ; The pressure drop across 
floats, as seen in figure 7 was nearly constant, except for the 
entrance region. The pressure drop due to float Ho.l was 
lesser than float Eos. 2 and 10. 

k) Effect of f loat-Height ; The heights tested in this study 
showed rotation-free performance in the limits of 
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D Float 
° Float 1 
6 Float 2 

4 - Float 3 



i o~r r r r i 1 r r r i r~ 

2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 


Flow rate in litre s/tninute 
FIGURE 7: TOTAL PRESSURE DROP ACROSS THE FLOATS 
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I n I o ic. 1 reeic; ing ox float = 0.8 001,1 
O for 1° taper rod (tod 1) 



Flow rate in li ter/ minute 


CALIBRATION FOR FLOAT 1 
FIGURE 8s EFFECT OF TAPER ANGLE OF ROD 
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4.2 Quantitative Results ; 

Multiple linear regression analysis was done to correlate 
tiie different dimensionless groups obtained through dimensional 
analysis. The form of the desired regression equation t-las 

Y i = B (°) + il boo ^(D + E i 0) 

i=i 

where p = no . of independent variables 

Ej_ = errors ■ 

B f s = unknown regression coefficients 

Three different cases were tried 

i) Regression among the logarithms of the given variables 

ii) Regression among the given variables 

iii) Regression among the logarithm of the depSndent variable 

and the given independent variables. 

Results were also obtained by deleting the independent variables 

which had very small regression coefficients. 

The best among all of the above was selected as the final 

regression equation. The selection was done on the basis of 

( 12 ) 

the following criteria 

i) The variance (RSS/Df) should be least. Here RSS is the 

residual sum of .squares, and DR is the degree of freedom. 

ii) The coefficient of correlation (R) should be the largest, 
o 

R indicates the percentage of variation in Y that can be 
accounted for by linear regression on the given independent 


variables 
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iii) The calculated value f ! , should be large 

, _ variance from li n ear regress lor ■ 

variance of residuals from regression 

When the value of I" is greater than the value obtained 
from the 3P- tables for the given degrees of freedom and 95$ 
significance level, then the hypothesis that no regression 
among the given variables exists is rejected. 

iv) for B(I)’s to be significant, the calculated. T factor 

for the respective B(I) should be greater than the value given 
in the T-tables at 95$ significance level, for the particular 
degree of freedom. 

v) The calculated Durbin-YYat son-statistics (D) should lie 
between the value DU and (4-DU) obtained from the D-tables 
for the specific degree of freedom. Then only the residuals 
will be independent of each other. 

The Durbin-Watson test can be carried out only for cases 
which at the maximum, have 5 independent variables, because DU 

values beyond that are not available even in the original 

( 19 ) 

paper; 

a) Correlation Among the Different Dimensionless Groups 
Obtained Through Dimensional -Analysis s In this case 


Y = (D2 - D3)/D 


XI = 

X3 = Hi/D 


Pw 


J WS - G 


"s 


8 


W' 


12 = H3/D 

X4 = D1 /D 
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15 = TMQ 
X7 = Cl /D 


X 6 = U 


A 


// w WS fft J> s - 


TABLE IV s DIFFERENT TEST-VALUES 

Number of independent variables = 7 

Total number of variables = 8 

Number of observation = 279 

Degree of Ereedoo = 271 



Log Eit 

Ordinary Eit 

Semi-Log Eit 

E ! 

1962.97 

1039.13 

833.87 

R 

0.9902 

0.9819 

0.9776 J 

Variance 

0.0029 

0.0002 

0.0066 

RSS 

0.7818 

0.0495 

1.7954 .I 


E* is much greater than the value obtained from the E-Tables . 
corresponding to a degree of freedom 271 . Hence -linear regression; 
in all the cases exists. 

R in all the cases is quite large, showing 96 - 98 $ linear variation! 
among the given variables. ; 

The variance in the case of ordinary fit is the least, 
hence it is the best among the different fits tested. 

Table V gives the coefficients of linear regression for 


the ordinary fit 
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TABLE Vs COEFFICIENTS OF BEQBESSION 
Degree of freedom for T = 273 

B(O) =-0.0983 


I B(I) 

Standard error 
in B(I) 

T-Factor 

1 1.0250 

0.0176 

58.1077 

2 O .2999 

0.0721 

4.1599 

3 0.0024 

0.0210 

0.1161 

4 0.0818 

0.0769 

1 .0628 

5 -0.4383 

. 0.3879. 

-1.1299 

6 0.5715 

0.1693 

3.3768 

7 -0.0411 

0.0197 

-2.0885 

From the students ’ t’ te 

st, for a level of 

significance 


0.2 

0.1 

Significant variables 

XI ,X2,X5,X6,X7 

XI , X2 , X 6 5 X7 

In sign if icant variables 

X3,X4 

X3,X4,X5 


Deleting the variables X3 and X4, and performing the regression 
analysis it was found that the value of R. decreased which 
implied that the new regression was no better than the previous 
one. 

The final correlation among the variables is 
| Y = -0.0983 + 1.0250 XI + 0.2999 X2 + 0.0024X3 + 0.0818X4 
I -0.4383X5 + 0.5715X6 - 0.0411 X7 (2) 




28 


"k) Corre lation for Pipe-Support Lengths In this case 



where M’ = slope of bottom calibration curve 
M = slope of remaining curve 
II = H4/D, X2 = M9, and 13 = Cl/D. 

Table VI depicts the different tests that were done. 


TABLE Yli DIFFERENT TEST-VALUES 



Humber of 

independent variables 

= 3 


Total number of variables 

= 4 


Humber of 

observations 

= 11 


Degree of 

freedom 

= 7 


Log-Fit’ 

Ordinary Fit 

Semi-Log Fit 

E ’ 

42.0833 

9.2501 

55.9146 

R 

0.9734 

0.8936 

0.9798 

Variance0.0322 

2.6089 

0.0245 

RSS 

0.2252 

18.2620 

0.1717 

D 

1.9089 

2.0686 

1.5956 


(13 1 ) 

The value of HU obtained from D-Tables v for a sample 


size of 11, independent variables 3, and an upper tail of 0.01, 
is 1.46. It is evident that in all the above cases the value 
of D lies between DU and ( 4 -DU) showing the residuals are 
independent of each other. Furthermore, the semi-log fit has 
the minimum variance alongwith a maximum value of R and F', hence 
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it is 

the best fit among 

the fits that were tried. 


The coefficients of 

linear regression for this 

fit are 

given 

In Table VII. 




TABID VII . C 

OEFEICIHIS OP.. EE GEES. SI Oil 



Degree of 

freedom for T =11 




B(0) = 5.4174 


i 

B(I) 

Standard error 
in B(I) 

I -Facto r 

i 

-0.6239 

0.0868 

-7*1459 

2 

-36.2462 

24,2510 

-1.4946 

3 

-7.2029 

1.3088 

-5 .5036 

Fro m 

the students 'T' test, for a level of significance 



0.2 

0.1 

Significant variables 

XI ,X2,X3 

XI,' X3 

Insignificant variables 

— 

X2 


To got a completely linear caliToration Y should have a value 1 t > 
then the final form of the correlation equation will he 

__ In (l)_= 30 + B1_ XI +22 X2 + B3 X3 (3) 

“ ~ ' — "i 

or 0.6239 (M/D) = 5.4174 - 7.202$ (Cl/D) ! 
i - 36.2 46? ! 


( 4 ) 



CHAPTER V 


Discussion 

linear calibration could not be attained In the lower 

end of the present flowmeter. This nonlinearity, was present 
in the initial 2 to 3 cms. of height. 


( 15 ) 


In the case of sudden contraction, Kaye and Rosen 


( 14 ) 


and 


Han showed, for polymer solutions and newtonian fluids 
respectively that, the pressure drop per unit length in the 
entrance region was comparatively larger than in the remaining 

portion. This pressure-drop had a linear relationship with tube 
length in both the regions. 

The present flowmeter, which is having a sudden expansion 
definitely has an entrance region where disturbed 
flow exists. The case here is contrary to that dealt by Kaye and 
Boson and W’5). Therefore , the pressure drop per unit 
length in the entrance region should be comparatively smaller 
than that in the remaining portion of the meter. The length, 

as the pressure drop profile, in this region will depend 

on the resultant , , 

rbarues caused by the sudden expansion 

and the presence of the tapered rod. 

The height of the float depends on the Impact heaa (l6 >, 

which in turn denprHo „„ 

on flow-rate. R 0 r a stationary float 

at equilibrium 

(P 2 -Pl ) Aa = V a g( j> a _ f ) (5) 

where p 2 = pressure at bottom of float 
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= pressure at the top of the float 
Ao = average area of cross-section of the float 
(P2~P-j) includes two effects viz. the. nearly constant pressure 
drop Ap', caused Toy the float and the pressure drop., AP, 
caused Toy the liquid flowing through the tapered annulus. 

The value of (pg-P-j)? and the calibration curve will depend 
mainly on the nature of variation in A p, because Ap' in its 
comparison, is quite small and constant. A greater impact head 
and hence a greater flow- rate will be required to lift the float 
to a particular height when At 1 is smaller, because in that 
case p,| , which is the pressure on float top, will be greater. 

The log-log behavior shov/n in an increasing order, by the 
calibration curve in the initial range (Pig. 5 ), is therefore, 
due to the entrance region which has a log-log pressure- 
drop— variation in the decreasing order. Hence,- to get a 
constantly linear calibration, it is necessary to avoid the 
entrance region, by varying the length of the metallic support. 

The effect of sudden expansion gets neutralized to some 
extent by the subsequent sudden contraction offered by the 
bottom of the tapered rod. Complete neutralization of the 
above effect can be achieved by adjusting appropriately the 
dimensions of the bottom diameter of the rod, and the pipe 
support. Smaller bottom diameters will provide lesser cons- 
triction and hence a lesser pressure drop. To have the same 
neutralizing effect, smaller bottom diameters will therefore 
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need longer pipe supports as compared to bigger ones, so that 
the pressure gets reduced to an amount that can he neutralized 
by the available constriction. It was seen that rods having 
bottom diameters more than 0.5 inches required very small 
pipe supports. 

The calibration line for a flowmeter having smaller taper- 
angled rod has greater slope than that having larger taper-angled 
rod of the, same bottom diameter. Hence, to bring the lower 
end of the calibration curve in line with the upper portion, 
it is necessary to provide an extra, pressure-drop than that 
required for n larger, taper angled rods. - ‘:s. A longer 
pipe support will do the needful. 

(17) 

For sharp-edged annular orifices, Bell and Bergelin 
showed that the discharge coefficient (C) , for a fixed Reynolds 
number decreases with an increase in the value of Z. Z is the 
ratio of orifice height, L, to annular opening, p=p. To get 
higher discharge coefficients it is necessary to have a low 
value of Z,. which can be done in the present case by either 

reducing or increasing — ^ — - . The restriction — jj — ss-0-3 

D2—C1 1 

limits the choice of — — . In a typical example where 

D2-A1 ^ 1)2 -Cl D2 r 

— j — = .5 , — |p— = . 01 , jp = . 6 and taper angle = 1 °, 

the height of the rod comes out to be about 15 cm. Therefore, 
D2-C1 

if |j — is taken much greater than . 01 , the range of the meter 

( 1 7 V 

will decrease appreciably. Bell and Berglins work showed 
that nearly for all Reynolds Numbers the orifice coefficient 



was minimum and remained constant for Z=0.1 to 0,4. Hence in the 
present case it will he desired to have SL^oi and ^,^2. 

Effective length as well as sensitivity of the meter depends 
on the tapper-angle of the rod. Both of these will he small for 
bigger tapper-angles. As explained hy Kaufman^®^ the taper- 
angle cannot he taken more than 8°, beyond which separation of 
flow occurs. In the given constraints the effective length of 
the rod for 1° taper will approxina tely he 15 cm. The sensi- 
tivity change in 1° ' taper and 0 o 5° taper is also considerable, 

Por taper angles more than 2° the effective length as well as 
the sensitivity will he quite lew* 

By increasing the float height it _is expected that rotation 
will start earlier than the limit of o,3 T However, greater 

float heights were not tested. 

It is believed that -the length of the vena-contraota, which 
depends on the annular orifice opening, is very small when 

— Q.3« Hence, the flowing fluid starts expanding in the 
float quite near its bottom. The smaller internal diameter or 
the greater float height disturbs the stream-lines near the 
surface, allowing them not to expand in their usual way. Because 
of this, it appears, some angular forces mqy come into action 
imparting rotation to the float. The exact reason of rotation 
is not clear. 



CHAPTER VI 


CONCLUSIONS M D RECOMMENDATIONS 

6.1 Conclusions s A flowmeter was constructed with a constant 
sided glass tube having a concentric tapered rod alongwith an 
annular float. The. float slided along the inner wall of the 
glass tube. 

Senstivity of the meter increased with a decrease in 
the taper angle of the rod. High senstivity can thus he achieved 
simply, Toy using a rod having a very small taper. In the 
conventional rotameter, this could have been achieved, only by 
resorting to the complicated design of Tarish!^ 

A linear calibration was found to exist between the float 
height and the flow rate. 

Heavier floats were more stable. 

The flowmeter showed a minimum of 99.5$ reproducibility. 
Rotation in floats was evident after (D2-D3)/D increased 
beyond 0.3. 

With regard to the design -of the flowmeter, it is concluded 
that the following constraints should he observed, 

i) (D2-A1 )/D <0.30 

ii) (D2-C1)/D ^ 0.02. 

iii) Cl/D 0.50 

iv) . DI/D ^ 0.78 

v) H3/D < 


0.01 
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vi) Hi/D <b_. 0*595 
vii) 0^2° 

By selecting the dimensions of the float within the limits 
of the above constraints,, .the weight of the float can be cal- 
culated. The flowmeter can then be designed for the given 
specifications by using the mathematical relationships obtained 
in Chapter IT. 

The fabrication cost of the flowmeter in the laboratory is 
estimated to be Rs. 25.0. This does not include the cost of the 
outer case and the permanent scale that will have to be used 
however, when this flowmeter is exploited commercially. The 
details of the cost estimation are given in Appendix. O . 

If this flowmeter is produced in bulk, the estimated 
cost of production of a 1-inch and a 2-inch diameter flowmeter 
will in any case not exceed Rs.75 and Rs.100 respectively. 

The present market rates of 1 inch and 2-inch rotameters are 
Rs.1200 and Rs.2000 respectively. 

The flowmeter is also capable of measuring opaque and 
translucent liquids. It is very cheap as compared to the 
present rotameters. 

6.2 Recommendations ; 

a) In Reference to flowmeter Design s It is recommended that 
the float, whenever possible, be made from materials having 
high density to ensure stability and rotation-free flow. Stainless 
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steel may be used as the material of construction. 

For very low flow-rates, glass tubes having appreciably 
small diameters, willjbe required, making it difficult to 
construct the flowmeter. This can be eluded by constructing 
floats out of thermosetting plastics, like; urea formaldehyde 
and cotton flock phono lies whose specific gravity lies in the 
range of 1 .3 to 1 . 5. 

The material of construction of the inner assembly has 
to be appropriately selected, from an economic standpoint, 
depending on the fluid to be metered. Except in the cases of 
concentrated acids and alkalis, thermosetting plastics may be 
used for the construction of the ipner assembly, other than the 
float. 

Using moulding processes to prepare floats and tapered 
rods, the necessity of calibrating individually each flowmeter 
of the same size and capacity can be eliminated. 

b) Recommendations for Future Study; The study of rotation 
in floats which is a general problem in all float-type flowmeters 
may be carried out. 

Experiments on the flowmeter using different liquids and 
glass tube diameters, may be made to check the validity of the 
mathematical relationships obtained among the different dimension 
less groups in the present work. 


*** 
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APEEUDIX 1 


TABLE . 1— As LET AILS OB BLOATS 


float 

No, 

Mat. of 
Constru- 
ction 

Wt.(WS) 

(g) 

<0 

•' s 
(g/cc) 

D1 V 
inch. . 

L2 ' 
inch. 

HI 

inch. 

H3 

inch. 

1 

A1 

5.7975 

2.715 

0.78 • 

0.605 

0.495 

0.040 

2 

A1 

6.9695 

2.715 

0.78 

0.600 

0.595 

0,054 

3 

A1 

6.0315 

2.715 

0.78 

0.604 

0.505 

0.075 

4 

Nylon 

3.0160 

1 .430 

0.78 

0.605 

0.490 

0.064 

5 

Teflon 

4.5680 

2.152 

0.72' 

0.605 

0.495 

0.069 

6 

A1 

7.1855 

2.715 

0.78 

0.605 

0.496 

0.054 

7 

A1 

7.0010 

2.715 

0.78 

0.601 

0.595 

0.065 

8 

A1 

5.7150 

2.715 . 

0.78 

0.631 

0.495 

0.061 

9 

A1 

5.6335 

2.715 

0.77 

0. 606 

0.480 

0.075 

10 

A1 

5.6235 

2.715 

0.78 

0.675 

0.495 

0.062 

11 

A1 

5.6940 

2.715 

0.78 

0.602 

. 0.495 

0.059 

12 

A1 

6.6805 

2.715 

0.72 

0.605 

0.495 

0.054 

Bloat 

Nos. 11 and 

12 had ' 

top taper 

angle of 30°, while in 

the 


rest of the floats it was 15°. 
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APPENDIX 2 

Values of the slope of the upper and lower portions 
of the calibration curves ior Eloat No«1, using different 


tapered 

rods and different heights of 

pipe supports 

are 

given in 

Table 2 .A 





TABLE 2 # A« SLOPES 

OBTAINED 

EOR DIEEEPENT 

PIPE 


SUPPORTS 




Values of Ml 

for Different H4 - 

M 

Rod No „ 

H4= 1 . 1 

H4=3. 6 

H4=4 . 7 

lv'X 

1 

0.32 50 

0.2060 

0.1440 

0.1000 

2 

0.3300 

0.1900 

0.1700 

0.0700. 

3 

0.4810 

0.2480 

0.1445 

0.0407 

4 

0.1700 

0.1489 

- 

0.1489 

7 

0.1530 

0.0715 

- 

0.0715 

- 

H4 = height of 

pipe support in cm. 



Ml = slope of the lower curve 
M = slope of the upper curve 
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APPENDIX 3 

MEIER CALIBRATION POP THE CASES OP DIFFERENT PLOATS AND RODS 

Initial reading of float for. all cases = 0.8 cm. 

Height of pipe support for all oases = 1.1 cm, 

HI = Height of Ploat in cm. 

Q = Volumetric flow rate in Litre/min. 


TABLE 3 -A; CALIBRATION EOR ROD 1 


HI 

Ploat 1 

Ploat 2 

Ploat 3 

Ploat 4 

Ploat 5 

Ploat 6 

Ploat 7 

Ploat 8. 

Q 

Q 

Q 

Q 

Q 

Q 

Q, 

Q 

1.0 

1.4458 

1.5333 

1.3932 

0.5934 

1 .0838 

1 . 7085 

1.5333 

1 . 6443 

2.0 

1 .7961 

1.8720 

1.6968 ' 

0.8678 

1.3348 

1.1172 

1.8603 

2.1581 

3.0 

2.0238 

2 . 1 1 72 

1.9304 

0.9262 

1 . 5330 

1.3741 

1.1055 

2.4675 

4.0 

2.1405 

2.2573 

2.0471 

0.9612 

1.6092 

2.5025 

2.2456 

2.6426 

5.0 

2.2456 

2.3741 

2.1581 

1 .0079 

1.7027 

2.6193 

2.3799 

2.8295 

6.0 

2.3390 

2.4792 

2.2573 

1 .0488 

1 .772 7 

2.7069 

2.4792 

2.8645 

7.0 

2.4441 

2.5901 

2 „ 3566 

1 .1247 

1.8720 

2.8061 

2.5959 

2.9930 

o 

CO 

2.5259 

2.7244 

2.4558 

1 .1947 

1.9420 

2.8820 

2.7302 

3.0513 

9.0 

2.6368 

2.8295 

2.5726 

1 .2473 

2.0238 

2.9930 

2.8178 

3.2190 

10.0 

2.7244 

2.9346 

2 , 6660 

1 .2940 

2.0997 

3.0455 

2.9171 


11.0 

2.8178 

3.0397 

2.7477 

3.3407 

2.1756 

3.1480 

3.0280 


12.0 

2.9229 

3.1480 

2.8412 

. — 

2.2456 

3.2440 

3.1480 


13.0 

3.0046 

3.2440 

2.9171 

1 .4341 

2.3099 


3.2190 


14.0 

3.1480 


3.0455 

1.5042 

2.4325 
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TABLE 3-B % CALIBRATION FOR ROD 2 


HT 

Float 1 

Q 

HI 

Boat 2 

Q 

HI 

Float 3 

Q 

HI 

Float 4 

Q 

HT 

Float 5 

Q 

1.0 

1.6795 

1 .0 

1 .8311 

1.0 

1 # 6268 

1.5 

0.9904 

1.0 

1 .2940i 

2.0 

2.0880 

2‘.0 

2.2048 

2.0 

2.0238 

2.0 

1 .0896 

2.0 

1.6443 

3.0 

2.3741 

3.0 

2.5492 

3.1 

2.3215 

4.8 

1 .2414 

2.9 

1.8311 

4.2 

2.5025 

4.2 

2.6777 

4.1 

2.4208 

5*6 

1.2765 

4. 6 

1.9654 

5.6 

2.6018 

5.2 

2.7477 

5.4 

2.5142 

6.7 

1 .3115 

6.2 

2.0296 

6.7 

2 . 6660 

6.1 

2.8120 

6.1 

2.5667 

8.5 

1.3699 

7.1 

2.0763 

7.2 

2.7010 

7.2 

2.8937 

7.1 

2.6426 

9.9 

1.4283 

8.3 

'2.1347 

8.4 

2.7769 

8.1 

2.9754 

8 . 6 

2.7361 

10.9 

1.4808 

9.8 

2.2048 

9.5 

2.8645 

8.7 

3.0280 

9.2 

2.7711 

12.0 

1.5509 

10.7 

2.2573 

11.5 

2.9930 

9.7 

3.1000 

11.2 

2.8995 

12.8 

1.6034 

11.7 

2.3040 







13.9 

1.6326 

12.5 

2.3507 


13.8 2.4208 
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TABLE 3-0 s CALIBEATIOH BOR ROE 3 


HI 

Bloat 1 

Q 

HI 

Float 2 

Q 

HT 

Float 3 

Q 

HI 

Float 4 

Q 

HI 

Float 5 

Q 

1.0 

1.8369 

1.0 

1.9829 

1.0 

1.7786 

2.2 

1.2531 

1.2 

1.4983 

2.0 

2.2865 

2.0 

2.4266 

2.0 

2.2398 

5.9 

1.3757 

2.0 

1.8253 

3.0 

2.6426 

3.0 

2.8295 

3.2 

2.5521 

6.7 

1.2816 

3.3 

2.1522 

4.2 

2.8120 

4.2 

2.9871 

4.4 

2.7244 

8.5 

1.4166 

4.4 

2.1989 

5.1 

2.8528 

4.9 

3.0397 

6.3 

2.8120 

8.9 

1.4341 

.5.4 

2.2223 

6.9 

2.9229 

5. 6 

3.0513 

7.3 

2.8528 

9.3 

1 .4574 

6.1 

2.2456 

7.6 

2.9462 

6.3 

3.1000 

8.3 

2.8879 

11.0 

1.5333 

8.2 

2.3157 

8.6 

2 ,9696 

7.8 

3.1715 

9.2 

2.9171 

12.2 

1 .5917 

10.1 

2.3624 

9.3 

3.0105 

8.2 

3.1960 

9.9 

2.9521 



11.1 

2.4033 

0.1 

3.0280 

8.8 

3.2441 

10.8 

.3.0046 



12.6 

2.4733 

2.1 

3.1480 



11.4 

3.0105 



13.8 

2.5200 


TABLE 3-D: CALIBRATION FOR ROD 4 


gloat 1 Bloat 2 Bloat 4 Bloat 5 Bloat 3 


HT 

Q 

HT 

Q 

HT 

Q 

HT 

Q 

HT 

. Q 

1.0 

0.2640 

1.0 

0.2090 

1.0 

0.1700 

1 .0 

0.1850 

1.0 

0.2325 

2.0 

0.4241 

2.0 

0.3949 

2.2 

0.2170 

2.0 

0.3073 

2.0 

0.3890 

3.0 

0.5759 

3.0 

0.5992 

3.3 

0.3014 

3.0 

0.4241 

3.0 

0.5525 

4.0 

0.7744 

4.0 

0.7744 

4.4 

0.3540 

4.0 

0.5408 

4.0 

0.7043 

5.0 

0.9378 

5.0 

0.9845 

5*4 

0.4124 

5.1 

0.6693 

5.0 

0.8794 

6.0 

1.1071 

6. 1 

1.7772: 

6*9 

0.5175 

6.0 

0.7919 

7.0 

1.2064 

7.1 

1.2706 

8.0 

1.4866 

8.7 

0.6517 

7.0 

0.9145 

8.2 

1 .3757 

8.1 

1 .4049 

9.3 

1.6910 

10.7 

0.8094 

9.1 

1.1655 

10.1 

1.6443 

9.2 

1.5625 

10.6 

1.8895 

12.6 

0.9437 

10.1 

1.2881 

11.1 

1.7727 

10.5 

1.7552 

11.2 

1 .9771 

13.8 

1.0196 

11.1 

1 .4049 

12.0 

1 .9012 

12.0 

1.9712 

12.0 

2.1172 

12.4 

0.9320 

12.1 

1.4107 

12.7 

2.0004 

13.0 

2.0822 

13.1 

2.2340 



12.6 

1 . 5684 

13.2 

2.0471 

14.4 

2.2748 

14.0 

2.2982 



13.4 

1.6443 

14.0 

2.1581 
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TABLE 3-E t CALI BRAT 101 FOR ROD 5 


HT 

Float 1 

Q 

HT 

Float 2 

Q 

HT 

Float 3 

Q 

HT 

Float 4 

Q 

HT 

Float 5 

Q 

1.0 

1.1305 

1.0 

1.1655 

1.0 

1 .0429 

2.1 

0.6226 

. 1 .0 

0.6342 

2.0 

1.2531 

2.0 

1.3582 

2.0 

1.2531 

3.8 

0.6693 

2.1 

1 .0254 

3.0 

1.4633 

3.0 

1.5100 

3.0 

1.3699 

4.5 

0 . 7044 

3.3 

1.1247 

4.0 

1.5567 

3.9 

1.6268 

4.0 

1.4750 

5.5 

0.7685 

4.3 

1 .2064 

5.0 

1.6501 

4.9 

1.7260 

5.1 

1.5720 

6.2 

0.8094 

5.8 

1 .2531 

5.0 

1.7494 

6.5 

1.9070 

6.5 

1.7319 

7.6 

0.8678 

7.2 

1.4283 

7.0 

■1 .8545 

8.2 

2.0997 

8.0 

. 1 .8778 

9.5 

0.9845 

8.7 

1.5567 

3 .0 

2.0588 

9.5 

2.2456 

9.6 

2.0355 

10.8 

1 .0429 

9 -9 

1.6618 

).1 

2.1639 

11.0 

2.3858 

1 1 .0 

2.1756 

12.2 

1.1247 

11.7 

1.8253 

1 .0 

2.2456 

12,4 

2.5084 

12.1 

2.3974 

13.8 

1.2064 

12.4 

1.8837 

’• .0 

2.3390 

13.0 

2.5551 

13.1 

2.3624 

13.5 

1 .1889 

13.8 

2.0063 

>.0 

2.4208 

13.8 

2.6543 

14.0 

2 .4500 






TABLE 3 -B: CALIBRATION BOR RODS 6 & 7 


pop 6 ' ■ ~ 3^7111 

Bloat 1 Bloat 2 Bloat 1 Bloat 2 

£ % — . — J P JE Q H T . 0 HT Q. 


1.0 

0.3715 

1 .0 

0.3598 

2.0 

0.5000 

2.0 

0.5058 

3.0 

0.6926 

3.0 

0 . 6926 

4.0 

0.8794 

4.0 

0.8678 

5.0 

1 .0429 

5.0 

1 .0488 

6.0 

1.1714 

7.2 

1.3874 

7.0 

1.3057 

8.7 

1 .6151 

8,0 

1.4574 

11.4 

2.0063 

10.1 

1.7377 

13.6 

2.3219 

13.2 

2.1114 

1 5.8 

2.5492 

15.5 

2.3741 

17.2 

2.7302 

17.1 

2.5492 

18.5 

2.8762 

19.6 

2.8120 

20.2 

3.0760 

20.8 

2.9229 

21.5 

3.1960 

22.6 

3.1240 

19.3 

2.9462 


1 .0 

1 .1363 

1.0 

1 .2181 

2.0 

1.3348 

2.0 

1.3582 

3.0 

1 .4574 

3.0 

1.5158 

4.0 

1.5275 

4.0 

1 .6151 

5.0 

1 .6151 

5.0 

1.7202 

7.6 

1 .8720 

6.0 

1.8369 

9.1 

2.0238 

8.8 

2.1114 

11.7 

2.2398 

10.7 

2.2923 

14.0 

2.4325 

12.8 

2.4792 

16.5 

2.6368 

14.0 

2.5901 

18.3 

2,8003 

16.0 

2.7653 

CO 

t> 

0 

CM 

2.9930 

18.2 

2.9462 

2.2.0 

3.1000 

20.0 

3.1240 

23.5 

3.2190 

21 .3 

3.2440 
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APPENDIX 4 ' 


DIMENSIONAL ANALYSIS 


As a first step to do dimensional analysis of the 
variables of the flowmeter it was necessary to identify those 
variables that affect the drag force R, acting between the 
flowing fluid and the stationary float. It has been seen 
that R is affected by the velocity V, density _P & the viscosity U 
of the fluid, and various geometric dimensions pertaining to 
the configuration of the tube and float, and to the orientation 
and position of the float. 

Among these dimensions a characteristic length, D, 
designating the diameter of the meter tube was selected. Using 
this 'characteristic length all other lengths were expressed 
as dimenionless ratios, XI, X2, X3, Xn,. 

Prom a practical stand point, it is easier to use the 
flowrate Q through the tube rather than the fluid velocity. 

Hence the former was used in the present analysis. The dimen- 
sionless length ratios are; 


Y1 _ D2-D5 7? _ 21 

A 1 - TQ 5 — D 


T7 _ Hi Y , _ is 

? X3 - -pj ? X4 — 


D ’ " D > x 5 - 2H2 


auAi , z6 = §1 


where D2-D3/D is a function of float height . 

Subject to the limitation of incompressible fluid, the inter- 
dependence of all the variables can be expressed in the following 
forms 
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F(R,Q, f , U, D, XI, X2,X3,X4,X5,X6) =0 (4.1) 

Equation 4.1 consists of 5 dimensional quantities and 6 dimen- 
sionless length ratios, expressible in terms of three basic 
dimensions, viz. mass length and time. In accordance with 
Buckingham's' ; Pi theorem, eouation 4.1 was expressed in 
terms of (5+6)-3 = 8 independent dimensionless groups. Equation 
4.1 can hence be written as 

f (Til, !7 2, XI, X2, X3, X4,X5, X6) = 0 (4.2) 

where 1 and 2 are two dimensionless groups consisting of 
R,Q, P , U and D, By arranging the above quantities in the 
form of quotient such that the physical dimensions cancel the 


following was obtained. 


TT i 


Q iJl 

D J R 


fi2 = TJ///R 

subsequent substitution of these values in equation 4.2 gave 

, XI, X2 , X3, X4, X5, X6) = 0 (4.3) 

If the consideration however, be restricted to a single float, 
at one given float position in a particular metering tube 



equation 4.3 ?/ill simplify to 


IT a TT 


( 44 .. ) 


(4.4) 


As stationary float is unrestrained in the vertical 


direction, the drag force must exactly balance the force of 
gravity. 
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Hence R = g 7 f ( f f - f ) (4.5) 

where g = acceleration due to gravity 
Vf = volume of float 
f = density of float 
= density of fluid 

Or, R = g m f (i' f -./’)/ f' f (4.6) 

here = mass of float. 

low, the unknown variable R of equation$«3) can be 
eliminated with the help of equation g-6) to give a relationship 
among the known variables, only. 




52 


APPENDIX 5 

MULTIPLE LIUBAB BESSES SIOH 

Basic Assumptions Underlying the_Estiaation Procedure s 

1. The expected valve of Y^, given X^, is a linear 
(in the parameters) function. 

2. The values of X selected for experimentation are nob 
random variables. 

3. The variance of residuals, equals the variance of 
Yb , and may be constant or vary with X. 

4. The observations of Y are mutually independent, i.e. 
the errors arc statistically independent. 

Based on those assumptions only the method of least 
squares yields unbiased estimators of regression coefficients B's, 
these will then have the smallest variances among the group of 
all unbiased linear estimators. It is therefore, necessary 
to chock the extent to which the above assumptions are valid. 

Here Yj_ is the dependent variable, while X’s are independent 
variables. 

Computational Procedure Adopted for Regression Analysis; 

Step-1 Preliminary Computations 
1. Selection and Transgenerations 

Subsamples are selected according to the specification 
on the selection card and data are transgenerated according 
to the codes. 
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2 . Sums : 

I 

sx(j) = y~~ i(i, j) j = i, — -Mi 

1=1 

where Ml = Total number ol variables 

H = Total number of observations 

3. Means ; 

I ’ 

X(J) = ( >— X(I,J)Xl, J = 1, Ml 

i=T 

4. Cross-Product Sums; 

__E 

S(K,J) = > X(I,K) X(I,J), J,K = 1,- Ml 

I=T 

5. Cross-Product of Deviations; 

J " . - " 

G-(K,J) = (X(l,K)-X(K)(X(l,J)-X(J) 

i'= i ... 

J ,K = 1 , Ml 

or Ct(E,J) = S(K,J) ~ (SX(K)*S-X(J)/ET 

6. Simple Correlation Coefficients; 

SC0SM(I,J) = r I-,J = 1, Ml 

(''12') 

Step-2 Regr ession Co mputations v y 

1. Inversion; 

C(I,J) = p x p Matrix, and is the inverse matrix of 
G(I,J) I,J = 1 p 

2. Regression Coefficients; 

JP ‘ 

b(i) = c(i,j) g(j,ii), i = i, — p 

j=i 
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■ If YG(J) = G(J,Ml) where J = 1, Ml. 

P 

Then B(I) = G(I,J) YG(J) 

J=T 

3. Intercept; 

P 

B(0) = X(M1 ) - > (B(I) XTIJ 

1=1 

Sum of Squares Atributable to regression 
P 

D = )> B(I) YG-(I) 

i=i 

Tests to check that the assumptions for linear 
regression do not get voilated. 

5. Sum of squares of deviation from regression; 

RSS = YG(M1 ) - D 

6. Coefficient of determination and multiple correlation 
Coefficient; 

R 2 = D/YG(M1) 

R =Jp? 

7. Variance and standard error of estimate; 

VEROR = RSS/ (N-p-1 ) 

S = J VEROR 

here (E-p-1 ) = Number of degrees of freedom = PE 

8. Standard deviation of regression coefficients; 

VCOVM(I) =./ VEROR • C(I,,I) 1=1, P 

9. T-values (to test the hypothesis that B(I)’s are zero); 

T( I) = B(I) /VCOVM(I) , I = 1, p 



10. 3? -Values (to test the hypothesis that there exists no 
regression) s 

I? =(D/p)/(RSS/PF) 

Step 5 s Test of Autocorrelation of Residuals v ; 

1. Error in calculated value of Ys 

"O’ (I) = Y-(B0 + B(l) ?Cl)+ — B(p)X(p)) 

I = i, .jff 

2. Residual sum of squares 

I 2 

UISQS = . U(I) 

1=1 

3. Difference between two consecutive errors; 

YDIFR = TJ(I) - n(l-l) 

4. sum of squares of the consecutive errors; 

_I_. • - 2 

DIEQS = > (U(I) - U(I-1 ) ) 

1=2 

5. Durbin-Watson - D-Statistics ; 


D = DIPQS/UISQS 
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AP PENDIX 6 

COST E SITUATION OF ONE FLOWMETER THAT 
WAS F A BRICATED II THE LAB . 

Total cost of Aluminium ref, trass pipes 


and brass plates PlS.5.00 

Cost of pyrex glass tube of 1 inch size Its. 4.00 

Cost of rubber corks Rs.2.00 

Approximate labour cost "Rs. 16.00 

Total i -Es. 21 .00 


Present rates 

Alu min iu m rods Rs . 9 . 0 0/Kg . 

Brass pipe and 

rod PLS.20.00/Kg. 

Pyrex glass tube 

(1 inch, size) Rs. 4* 00/Kg. 
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AP PENDIX 7 

C ALIBRATION OP PEE ROTAMETER 
The rotameter which was used as a means to indicate 


the flovi/ rate in the present flowmeter experiments was re- 
calibrated before use. The observations of the rotameter 
calibration showed a linear behavior. A straight line was 
fit to the data with the least squares method. Tests were 
made for 95 i° confidence limits to delete the outliers if any, 
from the least squares computations. However, no outlier 
was found to exist. The experimental data alongwith the 


computer program is given below on the next page. 

The final form of the least squares fit was 
HT = 8.5590(Q) - 0.6224 


Variance = 0.0306 

where Q = liquid flowrate in liters/min. 


HT = reading of float in cm. 


Variance = gealau al . & 9 U_ares 
Degree of freedom 


IiT 

Q 

1.9 

0.3040 

4.0 

0.5764 

8.0 

0.9967 

10.0 

1 .2192 

14.0 

1 .6783 

18.0 

2.1497 

22.0 

2.6519 

25.0 

3.0243 

25.3 

3.0490 

18.0 

2.1470 


HT 

Q 

2.9 

0.4528 

6.0 

0.7819 

9.0 • 

1 .1141 

12.0 

1.4564 

16.0 

1.9113 

20;0 

2.3949 

2 4..0 

2 ’.8841 

25.2 

3.0270 

25.5 

3.0713 


0 o o 
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H H 

i 

-I 


12 


10 

HT, cm j 


6 'i 


4 — i 


j.n 


O 

A 

X 


itial reading of float = 0.0 era. 

for pipe length = 1.1 cm. 
for pipe length = 3.6 cm. 
for pipe- length = 4.7 cm. 


ROD 2 
RIO AT 1 


d 

h 

* 


r 


P 


x 


.*/ 


<0 


<z 


/ 

/A 


/i 


„Jd(' 




■ 


-X 


0 


1.2 1.4 1.6 


; " i I 1 1 / 

1.8 2.0 2.2 2.4 2.6 2.8 3.0 


flow rate in liters/ minute 


FIGURE 9 s EPPECT OP PIPE SUPPORT LEITGTH 
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Initial reading of float = 0.0 cm. 



FI OKI EE 10 ? EFFECT Of PIPE SUPPOET LENGTH 
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14- 

12 ~ 

10- 

HT, off. 

3 - 

i 

6 - 

4 - 

2 ' 

0 — 

1.2 


Initial reading =0.0 cm. 


O 

A 

4* 



for pipe length = 1.1 cm. 
for pipe length = 3.6 cm. 
for pipe length = 4.7 cm. 


ROD 3 
PLOAT 1 



I 1 I !' 

•1 .6 1.8 2.0 2.0 



r i t '! i 

2.4 2.6 2.8 3.0 3.2 


Plow rate in liters/minute 


P1GTJRE 1 1 : EPPEC1 OP PIES SUPPORT LEIGTH 
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Initial reading of float =0.0 cm. 

20 — 




ooooooooooooonoooooor* 


THIS PROGRAM IS OF MULT 1 PLF-LINF AR-RFGRF $S ION 

a i =top nr a. of t * p c r r n n 1 inch. 

r ] =BOTTOM DIA . OF ROD, INCH. 

01= INSIDE 01 A. OF FLOAT, INCH. 

02= ORFICF nTfl. r-f FLOAT', INCH. 

03' =DIA.OF TAPER ROD C 0 R R F SPONGING TO HT , INCH. 
HI = LENGTH OF FLOAT, I’! r H. 

H2 = LFNGTII TAPERED ROE ,C'<. 

H3 1 ORFICF THICKM C SF, INCH. 

G= ECCL . DUF TG GR A V I T Y , C M / S E C • - S P C . 

FLCW = FLOW RATE, L/ W I ! :. 

HT = SCALF RFAOING OF FLOAT, CM. 

DEN$=D c: N$ I TY OF FLOAT, G/CC. 

DENW=n~HSI TY CF WAT"R , G v . /CC . 

MM = TOTAL NO. CF I N 0 . V A R I ARLES 
M] = TOTAL NO. CF VARIABLES 
M= DEGREE OF POLYNOMIAL 
H= TOTAL NO. OF nPS' : RVATIONS 
XY( J f Ml) • =DcP EMOFNT VARIABLE 
DIMENSION XY< 3) 

MM = 7 
J = 0 
G=981 . 

0=2.52 
DFNW=. 995 
U* .925 
DO 7 JM= 3 » 7 



READ (5, 501 )N5 

, NM) 

t C3 t 

A- 

■% 

4 , f 

H 2 

1 0 1 

FORMAT (21 

2»2F 

C ' ‘2: 

0 • -'J ! 

» F 6 . 

*} 

) 



DO I JM*J 

, NN 







REAL (5,2'’ 

0 ) D! 

»D 2 • 

t H I , 

w 

S , 

DEMSyH3 

200 

FORMAT (FA. 

-.2,2 

FF. 

? » F 6 

m 

A" £ 

H p 5 * ^ j 


N=N5 

R= (WS*G*(DEN$-DENW) )/DFNS 
FO=U/SQRT (OENW*R ) 

T A NA =2 • 5 4* { Cl — A I )/ (2.*H2) 

DO 1? 1=5. » N 
RF AD ( 5 ♦ 300 )HT , FLOW 
300 FORMAT (F5.1,F7.4) 

D3=Cl-{Ci“Ai) *CHT-.8) /-H2 
OOP =2.54* CD2-D3 I/O 

HOP ={ 10D.*FLCW)/(6.*D)*SQRT{Q£NW*DFNS/ ( W5*G*( DENS-DENW I) ) 
J= J+l 

*Vt Jr2t«2W v 54*H3/D 
XYf J, 30 =2. 54* HI, /O 
XY(J, 4)=?. 54*01/0 
XY(J,5)=TANA 
X Y ( J » 6 5 =F0 



xyu, ?)=;.. 54*ci/n 
xy< j,e )=oop 

12 COM THU r 

1 CONTINUE 
7 CONI T HUE 

INDEX*- 2 

00 4 1=1,3 

I NDEX~IMDEX+1 

CALL MURE !XY , J, MM, I N D F X > 

4 CONTINUE 
STOP 
■ END 

SUPROUT I Hp MLTRr ( XY , * , X * , I NDE-'X ) 

DIMENSION !40<% 5 ) ,XY( 400, «) , 5X( - Q ) ,G( F, P ),?(.?.), vr, (A) , VCOVM{ 8, 8) , 

1 SC CRM ( 8 , 8 ) , YF ST ( ~ 00 ) , T { . ■ ) , H ('•» GO ) , A S C OR M , P ) , C ( 8 , 8 ) 

TK = MM 

Mi=MM+l 

ITER=M1 

IF ( INDEX) 2?*, 5 55,2 '-4 
204 D0255J=1»MM 

002551=1., N 

2 55 X{ I , J ) =XY( I , J 5 

002561=1, N 

256 X { I , M3 3 = ALOG( XY { I ,M 1 ) ) 

WRITE! 6,607 ) 

607 FORMAT ( /2QX, 8 2HI NDr X POSITIVE , REGRESSION DF LOG Dl PENDENT VARIA8 
1LE0N GIV' ; N INCFPFNDJiNTVARI ARLES) 

GCTG565 

224 00260 J=l, Ml 

D0260I =1 , N 
XYMT =XY ( I J ) 

IF (XYMT >187,1 .88, 1.8? 

188 X C I , J ) =0. 

G0T0260 

187 X { I , J ) =ALOG f X YMT ) 

260 CONTINUE 

WRITE! 6,605) 

605 FORM AT ! / 2.0.X , 5 6H I NDEX NEGATIVE, REGRESSION ON LOGARITHMS OF ALL VAR I 
I ABIES/) 

G0TO565 

555 D0575 J= ! , M 1 

00575 1=1, N 
X ( I , J ) =XY ( I , J ) 

575 CONTINUE 

WRITF(6,606) 

606 FORMAT { / 20X ,4 1, HI NDEX ZERO, REGRESSION ON GIVEN VARIABLES) 

565 WR IT F! 6 , 0 511 -R,N,MM ■ 

80 FORMAT ( /I OX , I 6HN0. OF VARIABLES, I3/2 5X , 12H0B SERVATIONS 1 1 3/20X, 14HI 
1ND. VARIABLES ,12/) 

AA=0. 

0=0. 

, 0014 J= 1, Ml 

D015|=l, HI 

•' 

VCQVMf J,I )=0. 

■scormT:i,ui=o. 



15 COM T HUE 

B { J ) = 0 . 

S X ? J ) = ' . 

VG { J ) = ' . 

14 CONTI MUS 

TN = F LOAT { N ) 

DO 2'- J= l. Ml 
00201=1, M 

20 SX?J)=SX(J)+X(I»J) 

WR 110(6,900) 

9 A 0 FORMAT (/MX, 1 2HS.X = SUM OF INDIVIDUAL VAR . VALUES/) 

009C1 J = 1 , M 1 

90 1 WR I TF ( 6 » 19 ) J, SX { J ) 

19 FORMAT (/20X, 3HSX ( , 1 1 , ?! 1 ) = , F] 5.9) 

D011K=1,M1 
D01 1 J=K, Ml 
0011 1=1 

li G?K, J)*G(K, J) + XC I ,K )*X ( I, J) 

D0313K=1, Ml 
V D0313J=K , M 1 

313 G( K, J)=G(K , J)-{ ( 5X ( K ) *SX ( J) ) /TN ) 

D0I91 J=2,M1 
K*J-1 

D01 91 1=1 » K 

191 G{ Jj I )=Gf I , J) +0. 

WR I TE { 6,1-51 ) 

151 FORMAT ( /6X » 54HMA.TRI X OF SUM OF SQUARES AND CROSS PRODUCTS FROM MEANS/) 

IMS / } 

WRITE? 6,194)? (G ( T , J } » J = * ,MM), 1 = ] ,-MM ) 

194 F0RMAT{/6X,?F15. 6) 

001201=1, Ml 
0Ot2OJ*l,Hl 

SCORM ( I , J ) =G( I,.J)/SQRT(G(T,I)*Gf J,J) ) 

120 ASCORf I,J) = SCORM( I ,J ) 

WRITFU, 3? 2) 

122 FORMAT ( /29X »? 5HS ! MPLF CORRFLATION MATRIX/) 

WR IT F ( 6 , 1 24 ) ( { SC nRM ( I , J ) , J = l, Ml ) , I =1 ,M1 ) 

124 FORMAT ( /6X ,8F 1’.6) 

00551=1, N1 
YG(n = G(i, Ml) 

55 WRITE? 4 ,4 0)1, YG( I ) 

60 FORMAT ( /20X,3HYG ( , II ,2H)=,F15.6) 

CALL MATIN (G, NM, DETER ) 

WRITE ? A, 11 3 ) 

113 FORMAT? /3QX,1 5HINVERTFD MATRIX/) 

WRITFC6,! 14)? <G( I , J) , J=‘ , MM) , 1=1 , MM) 

114 FORMAT ( /6X ,TF 15.6 ) 

D042 J = 1 , MM 
00411=1, MM 

41 B( J)=SU) + G(J ,I)*YG!I) 

D=D + B? J)*YG(J ) 

42 CONTINUE 

R = SQRT f D/Y G ? H 1 ) ) 

RSS=YG?Ml)-0 
D’P=TN'“TK-1 * 1 
V€ROR=RSS/ OF 



WR I TF ( 6 , 64 ) YG ( Mj ) , V7R0R , DF 

64 FORMAT ( /4X , 1?HSQS. S'UMM D= , F 1 2 . 6 /6X , ?\ H RRDR=RF SI HUAI . S.S/OF=, 

1 F3 2.6/KX,9HnF FOR T=,F?Vi> 

WRITE! 6,61 ) 

61 FORMAT {/:?0X, HT , 5X, 4HR H ) , 5X, I6HST. ERROR OF B,S,4*,IHT) 

DQ6 ‘ T=’ , MM 
0066J=1,MH 

66 VCCVMU ,J) =GU» J)*VEROR 
AA=A A+B ( I ) *S: X ( I ) 

ORB= SORT ( V GOV M ( I , I }) 

T ( 1 5 = B ( I ) / ORR 

WR ITE{ 6,01 ) 1, 8! I) , ORB ,T ( I ) 

81 FORMAT Q7X , 14 tU, Fll • F, :X,F ! 5.^, IX, FI. X.F ) 

68 CONTINUE 

WRITE<6,71 } 

71 FORMAT ( /20X»? 6HVAR I-A.MCE COVARIANCE MATRIX/) 

WRITE! 6,72 ) ! ( VCOVM { I , J ) , .J = 1 , MM) , 1 = 1 , MM) 

72 FORMAT!/ if) X»? FI 5.6) 

AA = (SX(M1) — AA )/TN 
RSC=R*R 

F»(RSQ*DF) /{{ ] .-B5G)*TK) 

KM1=TK 

NMK=DF 

WRITE! 6 ,77 )D, RSS, AA , R , F , NMK 

77 FORM AT ( /"DX , ’ PHFXPL. SUM OF SOS. = , FI 5.6/20.X , 21HRFSI DUAL SUM OF SQS 
1 . = » F 1. 5 « 6/2 OX, 6 H B { 0 ) = , FI f>. 6/20X , ?HR= , FI 5 . 6/20X, 2HF= , FI 5.6/20X ,3HDF= 
1,14) 

WRITE16, *‘00 ) 

300 FORMAT <1H0,4 C X, 26HANA LYSIS OF VARIANCE TABLE //5X,1 9H SOURCE OF VAR I 
l AT I0N,.3X , 2HDF , 3 X , 1 ] HI UM OF SQS. » 3X» 1 2HC ALCUL ATED-F, 3X, 6HREMARK) 
WRITE 16,302 )K Ml , D, F 

302 EORMAT ( 1H0 , 13 X , 10HR2GP.S SS ION, 6X, J 3, F 11 . 6, 2X, F15 . 6) 
WRITE(6»304)N W K,RS$ 

304 FORMAT { 3 HO , 1? X, 8HRES T DUAL ,-6X » I 3, Fll . 6) 

NX =0-1 

WR I TFT 6 , *>0 6 ) N 1 , YG ! Ml ) 

206 FORMAT { 1H0 , 15 X, SHTOTAt , 6X, T3, FIX .6) 

WRITS (6,182) ■ 

18? FORMAT (1H0) 

YESTT=CU 

UISQS=0. 

DO 159 1 = 1, N 

no 159 J=1,MN 

YESTT = YFSTT+P ( J )*X( I , J ) 

150 CONTINUE 

YE STT=YE STT+A A 
XIMI=X! I, Ml) 

UT =X IMI-YE STT 
UI SQ=UI*UI 
UISOS=UISOS+U ISO 
YEST ( I. ) =YE STT 
U ( I ) =UI 
YESTT-O. 

159 CONTINUE 
01 FQS=0. 

• 00 '169 1 = 2,7) 



YD TFR=H ( I ) - U < I- 1 ) 

YD1 Fq=YPIFR*YCI. FR 
D I F 0 S = D I F 0 S+ Y C 1 F C 
! 69 CONTI MUf 

C n W O s T - r UJ R 3 I f W A T 5 n nj o- STATIST ICS 

DWDST=D I FO S/UI SOS 

WRITE ( 6 , 1.7 9 ) D tuDST » N » MK , UI SOS’ 

179 FORMAT! IMF , i; y, 27HDURRI ’• WATS HO D-$ TAT I FT I CS = , F 1 c . 6 /I 5X , I ? HOP. Si P.V-* 
1TI0NS = *I A /I D ,2' HF v Pl A O '- TORY V AR 1 A B !. r S = , 15/} 2X , f. HUT SQ$= , F? . h// } 

DO 1 8 1 J = ! » Ml 

no 1 R \ 1=1, N 

xn,j)=o.o 

181 CONTINUE 
C 988 RE ADA, MOOT A 
C I F ( NODTA ) 5 18, 183,1 

184 FORMAT ( 5X ,-FP . , F " . ' ,FF. ,F?. ■ , F- , F 7 . ' ,FP.' } 

RETURN 

END 

SUBROUT IMF MAT IN (A, N,D r TFR) 

. DIMENSION IPI VO( 5 ) ,A( 8, S) » INDEX! 8,2) , PIVOT (8) 

A B S F { X ) = A 5 S ( y ) 

C INITIAL I Z ATI CM 

10 DFTFR=1- 
M= 1 

15 DQ2C,J = 1 ,f| 

20 IPIV0!J)=0 
30 DOS 501=1 ,N 
C 

C SEARCH FOR PIVIT ELEMENT 

C 

40 AMAX=0. 

4 | 4$ ' I>01 05 J=1 » N 
50 I F ( I P IV 0 ( J ) - 1 ) 6 1,105,69 

60 DC3.00KM ,.N 

70 I F { I P IVO ( K ) -] ) 80, 100, 740 
PC- I F { A B S F { A M A X ) - A B 5 F { A { J , K ) ) 5 8 5 , 1 0 ' , X 0 0 
85 IRCW=J 
90 I C 0 L M = K 

95 AMAX = MJ,K) 

100 COM INU r 
105 CONTINUE 

110 T PIVr { ICOLM ) = IP I VGt ICDLMJ + l 
130 I F ( I ROW- 1 C OL M ) 5 4" , 2 6" , 1 < 0 
140 DETE R=-0 ' T FR 
150 D0?C0L=1 ,N 
160 SWA P=A ( I ROW , L 5 
170 A{ IROW»L)-A(I COL X, LI 
200 A { ICOLM , L ) =SW A p 
260 INDEX! I ,15 = I R C W 
270 I N C E X { I , = I C CL h ' 

31 > PIVOT ( T } = A { TC C L M * I COL K ) 

320 DSTER^pEff R*P IVOT ! I ) 

C 

C DIVIDE PIVOT ROW BY PIVOT ELEMENT 

C 

330 At ICOLM, ICOLM )=1. 



uu u u u u 


74' PC r--' 1=1, .J 

350 M ICOLM,L)=A( ICOLM f L)/PTVbT(I ) 

RE DUC.? MOM- PIVOT HOWS 

380 D055'~L1 = ' ,.N 
390 I F { L 1- 1 COL M H- CO , 5 5 0 , 00 
400 T=vfi (LI, ICOtM) 

4 20 M Li » ICflLM ) - : . 

430 DO 450 L=],N 
450 A (Li ,L) = M i: , L )-'• ( ICDLMt L)*T 
550 CONTINUE 

INTERCHANGE COLUMNS 

600 DO 710 I = I » N 
610 L = N + 1 — I 

620 IF{ INDEX (L »! ) -INDTXd. , 2 )) 630, 730,63 0 
630 JROW* INDEX (L, 1) 

640 JC0LM=INDFX(! ,?) 

650 DO 705 K=1 , N 
660 SWAP=A(K, JROW ) 

670 A { K , J R nw) = A ( K , J C 0 L M ) 

700 A(K,JCOLM]=SWAP 
705 CONTINUE 
710 COM I Mi IF 
740 CONTP'UF 

RETURN 
END 

■EEF-XIY 
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// JOB 

// FOR NAME 

♦NON PROCESS PROGRa v 

*IOCS< CAPD,TYP r WRI T EP. 3 

THIS PROGRAM IS OF ST. UMp FIT .BY LEAST SOS. 

95 PERCENT 0 CNF ID r: or F LIMIT IS TEST OH 
THIS PROGRAM WAS RIJM q*j jd.m J?00 

N = N 0 . OF OBSERVATIONS 
REAL M 

DIMENSION H{' C) , Q( 00) , 0 (20) , ERROR! 20) » GF.l ( 2 0 ) , Qf;2 ( 2 0 ) , K ( 2 0) , 
lPERORt 20) 

FO=£i. 

TANA=0. 

KQ-1 
K.P=2 
N= 1 9 

no i I = i j n 

READ (KP ,1 1 1 )H (I) ,Q(I) 

III FORMAT { F4. 1 ,1 X, F6.4) 

I CONTINUE 
11=1 

!F(II-1)8, 11*8 

8 KL=C) 

DO 9 I. = ), M 

DO 10 J=l, L 

I F ( I— K ( J ) ) 10* <3, in 

10 CONTINUE 
KL = K.L + 1 

H ( KL ) =H ( I ) 

O(KL) =Q(I) 

9 CONTINUE 
N=KL 

11 SUMH =0. 

SUMQ=0. 

SUNHH=0. 

SUMHQ=G. 

DO 2 1=1, N 
$UMH = SUHH+H( I ) 

SUMO=SUMO+Q ( ! ) 

SUMHQ*$UMHQ+H (I )*Q( I 5 
SUMHH=SUMHH+H ( I )*H( I) 

2, : S CONTINUE ■ 

PN* SUMHQ ) / ( SUMH*SUMH-PN*SUMHH 3 
C* I 1 . /PN ) * ( SU MO- M*SUMH ) 

SUME=0. 

DO 1 1=1, N 
Qci I )=M*H( T )+C 



- RR!,R { T ) =0 { i } -o' { T J 
SUK' = SUff'+ ! "RR.r& ( T ) 

PEROR ( 1 ) =c. RRf JR ( I )*; . / 0 ( I ) 

CONTINUE 

S JGMA=SQRT { SUM^/FLOAT 5 ) 

VE R = S UM P /F LQA T ( N- ] ) 

C4=C- 1.9^* SIGMA 
C2=C+] .9** SIGMA 

L =0 ‘ . 

DO 5 1=1, M 

QE J (I ) = M*H { I) +C« 

Qc2(.I )=M*H ( I ) + 02 
I F { Q { I ) -QE 2 ( I )) 31,-2 ,?> ■ 

31 I F ( Q ( I 5 -OF 1 ( I n'-,f ,5 
30 WR ITE (KO, 201) 1,0(1 ),H( I ) 

201 FORMAT ( 17, ?F? 2 . A ) 

L=L + 1 
K ( L ) = I, 

5 CONTINUE 
I F ( U 8 , 6 , 8 

6 WRITE ( KO, 30 1 ) 

301 FORMAT ( 5X , ?8H NO VALUE OUTSIDE 90PFRCEMT SIGMA LIMIT) 

W R I T E ( K 0 , 1 0 ■ ’ 5 C , M ,SJG M A , VE R 

100 FORMAT ( //8X , ?HC- , F8..A , "X » 2HM = , F 8 .4 , 5X, 5HSI GMA, F8.4, 
1 5X *3 HVFR , F 6 . &7 / ) 

DO 7 1 = 1, N 

WR I TE £ KO , 9 1 5H ( T> , Q H ) , Q ; ( I) , PEROR M ) 

91 FORMAT (AF1 5. ) 

7 CONTINUE 
CALL EXIT 

END 

// XEO NAME 

CCEMD 



